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Equations describing the propagation of radio waves in a horizontally stratified aniso- 
tropic ionosphere were developed by considering the limiting case of a large number of in- 
finiteshnally thin slabs of constant electron density and collision frequency. The quasi- 
longitudinal approximation was used. The pro|)agation equations appeared as four coupled 
first-order linear differential equations, coupled by gradients in electron density and collision 
frequency. The cjuasi-longitudinal approximation i)ermitted use of particularly simple 
forms for the coupling coefficients, these forms being amenable to simple analysis. Coupling 
between two ordinary or two extraordinary modes was found to be considerably stronger 
than cross coupling between ordinary and extraordinary modes. Cross cou})ling was rc^latc^d 
to the rate of change of the direction of the })hase normal. It was found that the reflection 
of VLF radio w^aves from the daytime ionos])h(;re is relatively insensitive to the angle of 
incidence on the ionosi)here except for highly oblique propagation. Whistler penetration 
was also found to be insensitive to the angle of incidence on the ionosphere. 

1. Introduction 

In this report, we shall consider the full-wave solutions for a very low frequency radio 
wave propagating obliquely into a horizontally stratified ionosphere. The effects produced by 
the presence of the geomagnetic field and by electron collisions will be considered in some detail. 
Although ionospheric reflection coefficients have been accurately calculated by Johler and Har- 
per [1962] for quite general conditions, the impetus for this study is to understand more fully 
the reflection processes of VLF radio waves and to understand the penetration of radio waves 
through the ionosphere. Wait [1960] and Barron and Budden [1959] have conducted extensive 
investigations of the reflecting properties of the ionosphere using the quasi-longitudinal ap- 
proximation and the mode theory; however, the actual reflection processes of the ionosphere 
have not been investigated. 

This report will consider the propagation of four separate magneto-ionic modes in a hori- 
zontally stratified ionosphere. The quasi-longitudinal approximation will be used in this 
analysis. The four modes have propagation constants given by the roots of the Booker quartic 
equation and represent two upward-propagating modes, ordinary and extraordinary, and two 
similar downward-propagating modes. 

The equations for the coupled modes will be developed through consideration of the iono- 
sphere as a series of step functions where, in between the steps, the electron density and collision 
frequency remain constant. Expressions will be written for the electric field intensity for each 
of the four modes within each step. These expressions will be related to the electric field in- 
tensity of a neighboring slab by means of the boundary conditions. Then the limit is taken 
as the width of each slab, and the change in electron density and collision frequency across each 
slab becomes infinitesimally small. Four differential equations are obtained for each mode, and 
they are coupled by gradients in the electron density and collision frequency. 

Although the general full-wave equations for ionospheric radio propagation are known 
[Budden 1961], they have not been studied in detail for oblique incidence because of the great 
complexity of the coupling coefficients. However, by using the quasi-longitudinal approxima- 
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tion, reTaarkable simplifications can be made which permit the simple analysis of the behavior of 
VLF radio waves in the ionosphere. Limitations on the quasi-longitudinal approximation will 
also be discussed. 

The reflection processes will be discussed as the coupling of the upgoing modes to the down- 
going modes. Explicit representations of this coupling will be given as a function of altitude 
for a typical daytime ionosphere. The reflectivity of the D region as a function of the angle of 
incidence will be qualitatively indicated. The penetration of very-low-frequency signals 
through the ionosphere will also be qualitatively discussed by examining how much of the wave 
is coupled into the downgoing modes and how much is attenuated when passing through the 
ionosphere. The optimum conditions for penetration from the ground into the ionosphere will 
be indicated. 

2. Development of the Propagation Equations 

2.1. Basic Considerations 

In this section we shall begin by showing the equations for a wave propagating in a uniform 
magneto-ionic medium. We shall then express the boundar3^ conditions for a wave incident 
upon a plane boundary where a discontinuous change in the propagation constant occurs. 

The quasi-longitudinal approximation for the propagation constant is given by 

^=1(^2+ ""^ — -Y (1) 

C\ — co+^^^±a)^ cos ^/ ^ 

where 

co = 27r/, 
cx)p= (square of the angular plasma frequency), 

J/ = electron collision frequency, 
o)m= — - electron gyrofrequency (co^ is negative because the electron has a negative charge), 

m = electron mass, 

6 = electron charge, 
iVe= ionospheric electron density, 
5o=magnetic induction of the earth's field, 

^= angle between direction of wave normal and direction of the geomagnetic field, 

c= velocity of light. 

Equation (1) is valid in the limit that [Budden, 1961, p. 119] 
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Suffice it to say for the present that (2) holds for VLF propagation provided that 6 is not ''too 
near'' 90 degrees. Equation (2) will be investigated in more detail later in this paper. 

Returning now to (1), we see that because of the presence of a magnetic field through the 
term co^, there are two modes of propagation and the medium is also anisotropic. Let us denote 
k of (1) as ki when the plus sign is used and fe when the minus sign is used. We notice that if 
the conditions of (2) are satisfied and if cos 6^0, a wave of propagation constant ki will gen- 
erally propagate and a wave of propagation constant k2 will be heavily damped. 

Before writing analytic expressions for the propagating waves, we shall adopt the coordinate 
system shown in figure 1. The magnetic field Bq is confined to the x-z plane and is inclined at 
an angle 8 to the a:;-axis. The electron density and collision frequency is to vary only in the z 

direction. The direction of propagation is denoted by the vector U at an angle \l/ to the 2:-axis 

and with an azimuth al angle </>. The unit vector 1^ is therefore given by 

-> ^ -^ -> ^ ^ 

1a;=12 cos \l/-\-(lx COS 0+lj/ sin 0) sin \p, (3) 
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Figure 1. Coordinate sy stein. 
The geomagnetic field vector may be represented by 

Bo=Bo{-lz sin 5+t^ cos 5) . (4) 

It is now easy to compute cos in (1). It is given by 

cos ^=— cos ^ sin 5+sin x// cos d cos (j). (5) 

In the quasi-longitudinal approximation, mode 1 will be right-handedly circularly polarized 

-> -> 

and mode 2 will be left-handedty circularly polarized. Let ly^ and Ig^ be two vectors forming 

-> 
a right-haiided orthogonal coordinate set of unit vectors with 1;^- Then, in general, the electric 

vector of a plane-propagating upgoing wave in a homogeneous medium may be represented by 

-^ -> -> 

E''(x,y, z)=Ei{ly'-\-il,^) exp {—iki [cos 1/^1^+ sin i/'i(?/ sin </)+x cos (/>)]} 

+E2(ly' — il2') exp {—ik2 [cos i/^s^+sin \l/2{y sin <^+x cos </))]} (6) 

where Ei and E2 are constants. 

-> 
Let us further define 1^/ as being parallel to the x?/-plane; then 



lyf = ly cos </)— Ix sin 0, 

-> 
so that 1,/ may be determined; namely, 

-^ -^ -> -^ -^ ^ -> 

l2, = l^,Xl2/' = l2 sin \p—(lx sin 0+1;^ cos 0) cos \(/. 



(7a) 



(7b) 



Therefore, when l/^ points along the x-axis, tlien 1^/ is along the y-axis and l,- is along the 2-axis 
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2.2. Reflection From a Plane Boundary 

When a wave is incident on a boundary between two regions with differing refractive 
indices, part of the wave will be transmitted across the boundary and part will be reflected. 
The reflected wave corresponds to two more modes with propagation constants. 



c L — co+^ 



CO COp 



iv± (j^ra (cos -^ siu 6+sin ^ COS h cos <^) 



]' 



(8) 



where we shall say Ats is the mode corresponding to the minus sign, and h^ is the mode corre- 
sponding to the plus sign. Equation (8) differs from (1) in that the sign on cos ^ has been 
changed. (Wherever cos -^ or sin ^ appears, it will always be understood that the real part 
is positive.) 

The electric field of the downward propagating waves is given by 

-> -> -^ 

E^=E^(\y'-\-i\z^ exp {%kz [cos i/^gS+sin 1^3 {y sin 0+a^ cos </>)]} 

-> -> 

-\-E^{\.y>—i\z') exp {%k^ [cos ;/'4S+sin )/'4 {y sin 0+x cos 0)]} . (8a) 

When using this expression, it must be remembered that k^^ is the propagating mode and kz 
is the evanescent mode, under the condition that the sign on the projection of the direction 
of propagation on the geomagnetic field, cos ^, has changed from that appearing in the upgoing 
modes. This will usually be the case for VLF propagation in temperate and polar latitudes. 
When cos d does not change sign between the incident and reflected modes, mode 4 is evanescent 
and mode 3 is propagating. 

Let us now consider a plane boundary in the x-y plane located at Zi. Let the propagation 
constant at z<C^Zi be k\^ and at z^Zi^ k\^^ where j = 1,2,3,4, corresponding to the two reflected 
(downgoing) and two incident or transmitted (upgoing) waves. Further consider another 
boundary at z—Zij^i such that both above and below the boundary at Zi there are two upgoing 
and two downgoing waves. 

The boundary conditions are that the tangential component of £'and of vX^'be continuous 

across the boundary. Therefore, by equating the x and y components of E and VX-B we 
obtain four equations relating E] to E]'^^. Thus the boundary conditions can be expressed in 
the form 
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(9) 



where 



and where 



Ej=Ej exp { —ik'j[e(j) cos T/^j^^+sin \[/'j{x cos 0+^ sin (j))]} , 



-^ -> 
The elements of the square matrices are computed from the x and?/ components of {ly'±ilz^) 

and from the x and y components 



hjX {ly ± ilz')kj^ {I,' =F ily'kj 
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The matrix elements, alj, may be written 

a}y=— siii(/)+i(— l)-^e(j) cos (j) cos 4/ J 
aij= cos (l)-{-i{—l)^e{j) sin cos ^j 
alj=[—€{j) cos cos \^i+i(— 1)^' sin 0]ty 
a4y=[ — e(j) sin </> cos \p'j—i{—l)^ cos <^]tj-. 



(11a) 
(lib) 
(lie) 
(lid) 



In order to satisfy the conditions of (9) at any point (x,y) in the z=Zi plane, the auxiliary 
condition 



k^j sm\l/j=— sini/'o 



(12) 



must be satisfied for all values of j and at all levels in the ionosphere where \l/o is the zenith 
angle of the wave below the ionosphere. Thus when an electromagnetic wave enters the 
ionosphere it will split into four modes, two transmitted and two reflected, all propagating in 
different directions. However, the azimuthal direction for the phase normals must be the 
same for all four modes. 

The problem of determining the angle ^^ is by no means trivial because from (1), (5), and 
(8) it can be seen that kj is an irrational function of sin xj/]. The solution of (12) for sin yp] or 
cos \l/] is equivalent to solving the Booker quartic equation for the four roots k] cos \p] given 

the value of — sin \po. However in practice, a simple iterative solution for sin \l/j can be found. 

For example, if sin \p^j is known and it is desired to find sin \//]'^^ then we let sin i^jy be the Zth 
step in the iterative process described by the equation 



siniAj.V- 



CO . 

— sni if/Q 
c 



(13) 



^j+'(siin^lV-,)' 



where 



siniA'/o'=sini/']. 



This simple routine will converge in a very few steps. 
Returning now to (11) we notice that 

—alj sin 0+a4j cos <^=1 

alj cos (t>+alj sin </>=(— l)"^'e(j) cos xp] 

alj sin (j>—a\j cos <j)= {—\)^k] 

alj cos (l>-\-aij sin 0= — e {j) cos i/'jA:j- 



(14a) 
(146) 
(14c) 
(14^) 



which results in a considerable simplification over (11), so that (9) may be written in the form 
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-kl + ' COS ^^ + ' fc|+' COS ^1+' 4+' COS ^l+\ 

where Ej is now understood to mean Ej exp [—ie{j)k} cos^js,]. 
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(15) 
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2.3. Transition to the Continuous Ionosphere 

In the preceding section, we considered plane waves incident on a sharp boundary be- 
tween regions of different electron density. Let us now examine a large number of steps 
which approximate continuous ionospheric conditions in a manner shown in figure 2. Let 
each of the intervals in figure 2 be of width A^ and let the index "i^' in (15) refer to the ^th 
interval. Then the quantities in the (^+l) interval may be represented in terms of the ith 
interval, provided that the change between the intervals is small, by the relations 



^ dz 
COS \^}+^=cos i^i+-T- (cos 1^5) A3 



(16a) 
(166) 
(16c) 



and 



:^xp [-e{j)ik)^' cos)Ai+'^z]-exp [-e{j)ik) cosiAj^,]X ^l-^6(i) ^ ik] cosi^D^.AsI- (17) 




ALTITUDE 



Figure 2. Representation of a continuous 
ionosphere as a step function. 



Upon combining (16a) and (17) and neglecting terms of order (A^)^ we get 

Ei^' = E) exp [-ie{j)ki cos ^J3,]+exp [-ie{j)k) cos r,z,] ^-ie{j)Eiz, ^ {k) cos \^})+J^ E]\ Az 

which may be written in the more desirable form, 

n+'=Ei4-\'S Ei 



E'j-^^=Ei+^-g Ei+u{j)E]k) cosv^ijAz. 



(18) 



Upon substituting (16b), (16c), and (18) into (15) and using the variable z instead of the 
subscripted variable Zi we get a coupled set of differential equations which express the bound- 
ary conditions between layers infinitesimal thickness Az\ 
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(19) 



This equation is somewhat awkward in that the derivatives of the Ej are mixed. The 
form of these equations can be considerably simphfied if we note that the first term on the 
right-hand side of (19) can be factored 
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(20) 



where Q is the diagonal matrix. If we let E be the colunni matrix containing Ej then (19) 
has the form 



AE' = {AQ-A')E 



(21) 



where the prime denotes differentiation with respect to z. The form of (21) may be consider- 
ably simplified if we multiply through on the left by A~^ obtaining 



E' = QE-A-'A'E. 



(22) 



Since Q of (22) is diagonal, we see that the elements of the matrix A~^A couple the modes 
together. We further see that coupling is dependent upon gradients in the electron density 
and collision frequency. If the coupling matrix is small, then all modes propagate almost 
independently and 



E 



r^exp [—ie^j) kj co^ \l/jdz]. 



(23) 



If we let (A ^)ij represent an element of A ^ then the {%, j)th element in the coupling matrix 
A~^A=C vcidij be written 

0,j=e{j^f^{cos^M-^)KA-')2i-h{A-'),.]-j^^^^ (24) 

Here we notice that the coupling between modes is produced by gradients in electron density 
and collision frequency, and also that when the ionosphere is homogeneous, the modes propa- 
gate independently with solution given by (23). 
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3. Coupling Matrix 

The coupling coefficients corresponding to the matrix C have been computed by Heading 
(see Budden [1961]) for the general case corresponding to the full Appleton-Hartree formula 
without approximation. His expressions for C are quite complicated and difficult to interpret 
except by computing numerical examples. The use of the quasi-longitudinal approximation 
permits considerable simplification of the coupling matrix to the point where physical inter- 
pretations are possible without making use of detailed computations. 

In order to simplify the notation used here, let us write 

C; = COS \pj. 

Then we may begin the computation of cy by finding the determinant of A, namely 

\A\=-{h^-h){h-^h){cr+Cz){Ci+Ci)-{h—kz){k2—h){c2-Ct)ici-c^. (25) 

A~^ is then computed by finding the cof actor of A, making use of (25). The elements of C 
are then found by substituting A"* into (24) and rearranging terms; thus 

\A\Cn=c[{{kx-]-ki){k2+kz){c2+Ci) — {ki-kz){k2-h){Ci—Cz)] 

+ ^1 [ (^2 + ^3) (ci + Cz) (C2 + C4) — (fe — ^-4) (ci — Ci) (c-i — C3) ] (26a) 

\A\Cu=C2{h-ki){h+k3){c3+c,)+k2ikz+ki)ic,-c-d(,C2+cd (26b) 

\A\Ciz=-C3ih+k4){k2+kz)ic2+Ci)+k3{k2-ki)ic,-C3)(c3+c,) (26c) 

\A\Cu=Ci{k3+k,)(k2-ki){c2-C3)-ki{h+kz){c2+Ci){c3+Ci) (26d) 

\A\C2i=c[{h-ki)(h+k,)liCi+Ci)+k'iik3+k,)iCi+C3)(ci-cd (26e) 

\A\C22=C2[{k2+k3){h+ki){ci+C3) — ik2-ki)iki-h)ici-Ci)] 

+k2[(h+k^)(c2+c,)(ci+C3)-{h-kz)ic2-C3){ci-Ci)] (26f) 

\A\C2i=c^{k3+k^)(h~■h)(Cl-c^)-k3{kl+ki){c3+c^)(Cl+c^) (26g) 

\A\C2i=-c[(kz+ki){ki+k,)(ct+C3)+k[{ki-k3){cz+Ci)ici-Ci) (26h) 

\A\C3i==-c[{ki+ki){ki+k2)(c2+c,)+k[{k2-k,){ci+C2){ci-c,) (26i) 

\A\Cz2=C2(k2-ki)ih+k2)(ci-Ci)-k2{ki+ki){ci+C2){c2+c,) (26j) 

\A\a,=4(k,+k,)(k2+h){c2+c,) + (k,-k,){k2-kMc:-Ci)] 

+ks[ih+k\){ci+C3)(c2+c,) + ik2-kdic2-C3)(ci-Ci)] (26k) 

\A\C3i=-c[{ki+ki){k2-ki)(ci+C2)-kiiki+k2)ic2+Ci){ci-c,) (261) 

\A\Cii=^c[(ki+h){k,--k3)(c2-C3)-k[{k2+k,)ic,+C2){ci+Ci) (26m) 

\A\C42=---C2{h+k2)ik2+h){ci+Cz)+k2(,k,-kz)ici+C2){Ci-Ci) (26n) 

\A\Ca=-c'3ik,-ki)ik2+k3){ci+C2)-k3(ki+k2)ici+C3)ic2~C3) (26o) 

\A\Cu=c'A{ki+hd{k2+kd(ci+Cz) + {h~k,){k2-ki)(c2-e,)] 

+k'i[ik2+k3)ic2+Ci)(,c,+C3) + iki~k3)(ci-c,) (02-03)]. (26p) 

The form of (26) is such that the effects of approximations and special cases can readily 
be seen. In the case of vertical propagation, we know that Ci=C2=C3=C4=l, c.=0, and 
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]ci=ki=kjc and ^•2=^3=^0- We tluis get the simple form for the coupling matrix 
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(27) 



From the above form of C we see that an upgoing ordinary wave is reflected into a downgoing 
ordinary wave and that an upgoing extraordinary wave is reflected into a downgoing extraor- 
dinary wave. We see that modes 1 and 3 propagate independently of modes 2 and 4. 

In order to discuss other special cases and approximations for oblique incidence, it must 
be kept firmly in mind that modes 1 and 4 are propagating and modes 2 and 3 are nonpropa- 
gating or evanescent. In the light of our previous discussion, this means that cos d in (5) 
must, upon reflection, undergo a change in sign. 

A simple example is propagation in a plane perpendicular to the magnetic meridian. 
Under these conditions cos (/> of (5) and (8) is zero, so that Ci=c^=Cj:, ki—ki=k^, c2~Cs=Cq, 
Sindlh=k^=ko. As a result of these simplifications, the coupling matrix becomes 
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(28) 



The form (28) will be identical for the case of a vertical magnetic field and oblique propa- 
gation in any direction. We notice that C2 has a form identical to that of (27) and that the 
effects of cross coupling between ordinary and extraordinary modes on oblique propagation in 
the plane perpendicular to the magnetic meridian are proportional to the rate of change of the 
cosine of the zenith angle of the phase normal. 
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In order to gain further insight into the properties of the coupHng and thus the interaction 
of the various modes of propagation in a nonhomogeneous ionosphere, we shall calculate the 
values of kj and Cj for a typical ionosphere. Figure 3 shows the electron density of a normal 
daytime ionosphere for middle latitudes along with a collision frequency. Daytime conditions 
are chosen so as not to impose too much of a violation on the quasi-longitudinal approximation. 
We shall also assume a constant magnetic field strength of 0.45 X10~* weber/m^, a dip angle of 
60 degrees, and a propagation frequency of 10 kc/s. The values of kj and Cj as a function of 
altitude are computed from (1), (5), (8), and (13) and the magnitudes of kj and Cj are shown in 
figure 4 {a-K) and the arguments in figure 5 (a-h) for several initial values of yp and </>. Figures 
4 and 5 show that the conditions imposed by propagation in a plane perpendicular to the mag- 
netic meridian represent generally valid approximations for propagation in any direction be- 
cause of the very small differences between Ci and C4, and C2 and C3. 




ALTITUDE, km 

Figure 3. Model ionosphere showing electron 
density and electron collision frequency. 

The second term in (28) is quite easy to interpret because it is made up of logarithmic 
derivatives of the propagation constant. Let us write k^^x in the form \kQ,xW'^ where \k\ is the 
magnitude of k and 7 is its argument, shown in figures 4 and 5, respectively. Then the nonzero 
matrix elements of C^ in (29) may be expressed as 

(7?,-l/2ln'|/:| + l/2i7o... (29) 

Thus 7' is the slope of the curve on figure 5 {a-h) and In' \k\ is the slope of the curve shown in 



figure 4 {a-h) since \k\ is plotted on a logarithmic scale. Similarly 
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Figure 4. Magnitude of the propagation constant and of cos \f/. 
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Figure 4. — Continued. 
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Figure 5. Arguments of the propagation constant and of cos \p. 
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nonzero elcinents in C^ of (28), may to a good approximation be written in tlie form 



Cq'I'Cx 



d/2ln'|c| + l/2i7;, 



m 



where yc is the argument of the cosines. Equation (30) is valid to within about 10 percent 
because the magnitudes and arguments of Cj are approximately equal. 

The coupling that produces reflection between the two ordinary modes and between the 
two extraordinary modes is the strongest and also it is relatively independent of the angle of 
incidence from below the ionosphere. The coupling between modes 2 and 3 is stronger than 
between modes 1 and 4 because modes 2 and 3 undergo a large change in argument between 
the altitudes of 65 and 80 km.^ 

Cross-coupling terms shown in matrix C^ are relatively weak except for highly oblique 

incidence. Factors of the type -^ — - have a magnitude that varies between 1 and I/V2 and 

an argument that varies between and 7r/4. Thus terms representing coupling between modes 
1 and 3, and modes 2 and 4 are characterized by the logarithmic derivative of the c^s. However, 
reflective coupling due to the elements of C^ occurs primarily between the altitude of 65 and 75 
km, whereas reflective coupling due to the elements of C^ is more in the 70 to 80 km range. 
Thus careful observation of the polarization of the reflected wave should allow observation of 
two different regions of the ionosphere. 

Coupling between the two upgoing modes and two downgoing modes, represented by the 
matrix elements {Clz, Cl^ and (6^34, C43), contain factors of the type {ko—kx)/2ko,x' The 
magnitude of this factor varies from zero at low altitudes, to 1/^2 at high altitudes. Thus it 
is to be expected that coupling between the upgoing modes or between the downgoing modes 
be somewhat less than reflective coupling between ordinary and extraordinary^ modes. 

4, Some Comments on the Conversion of Spherics Into Whistlers 

As a wave penetrates the D region of the ionosphere, some of the wave energy is converted 
into other modes, which results in tlie process of reflection, while some of the energy is dissipated 
by electron collisions with the neutral particles. In order to obtain an estimate of how much a 
penetrating wave is attenuated, we shall compute the attenuation coefRcients of the WKB 
solution of mode 1 given by 

Ei=Eo exp —i ki cos \pi dz \y 

where the attenuation coefficient is given by 

exp— lui{kiQ,o^\l/i)dz\ 

where Im {x) denotes the imaginary part of x. Table 1 shows some attenuation coefRcients of 
mode 1 for several initial values of ^ and (/>. 

Table 1. Attenuation of the WKB solutions for the upgoing propagating mode as a function of the angle of incidence 



Angles of incidence 


Amplitude 


Energy 






attenuation 


absorbed 


detj 




Nepers 


db 


<t> = \80 


^ = 15 


1.52 


13.2 


= 180 


=80 


1.54 


13.4 


= 180 


=45 


1.59 


13.8 


= 180 


=60 


1 71 


14.8 


= 180 


= 75 


1.91 


16. 5 


= 135 


= 15 


1.56 


13.5 


= 135 


= 30 


1.58 


13.7 


= 


= 15 


1.72 


14.9 



' Figure 5 (a-h) is scaled such that a distance of one radian is the distance represented by a factor of e on figure 4 (a-h) . 
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It can be seen from table 1 that the total absorption of, say a whistler signal, is relatively 
independent of the angle of incidence. It increases only gradually with increasing zenith 
angle. The favored azimuthal direction is in the same azimuth of the magnetic field. The 
optimum angle of incidence for minimum attenuation is somewhere between the zenith and 
the direction of the magnetic field. The small variation of total attenuation as a function 
of the zenith angle (much less than a sec xp variation) results from the sharp change in the 
direction of the wave normal toward the zenith as the refractive index of the ionosphere increases. 

Also, as indicated in the previous section, the conversion of mode 1 into the downgoing 
modes, 3 and 4, is somewhat less for more-nearly-vertical incidence. Therefore we may con- 
clude with some reservation that the optimum direction for the penetration of very low fre- 
quency energy into the ionosphere is that of a nearly vertically-propagating wave; however, 
the angles of optimum penetration have a rather broad maximum. 

The complete picture of the injection of whistler energy into the ionosphere must await 
the complete solution of the coupled equations shown in (22) . However, from figure 4 (a-h) 
we see that, as far as the direction of the phase normal is concerned and regardless of the 
initial direction of the incident spheric signal on the ionosphere, the whistler is injected into 
the ionosphere with its phase normal propagating vertically upward. Therefore, once the 
electron density gradient above 105 km flattens out, the ray path as shown by Storey [1953] 
will have a direction given by 

tan a^ — % cot d 

where a is the zenith angle of the group velocity vector. Therefore the ray path will assume 
a direction between the zenith and the magnetic field direction above the E layer. Figure 6 
shows some probable ray paths in the magnetic meridian generated by a lightning discharge. 
Since with a vertical dipole source most of the energy is radiated obliquely into the ionosphere 
and since reflection and absorption of mode 1 is not critically dependent upon the zenith angle 
of injection into the ionosphere, whistler sources may be quite distant from their apparent 
sources determined by tracing magnetic lines of force. Studies of whistler propagation [Helli- 
well and Morgan, 1959] have shown this to be the case. 




SPHERIC SOURCE 



Figure 6. Probable whistler ray paths through 
the ionosphere. 

5. Limitations on the Quasi-Longitudinal Approximation 

The great simplifications made in the analysis were possible because the very-low-fre- 
quency wave was assumed to be circularly polarized. If this were not the case, then the 
matrix A would not have had its simple form, and its inverse would have been even more 
complicated so that the coupling matrix A~^A' would have been very difficult to analyze. 

The condition that the quasi-longitudinal approximation be valid is that 

,^,o co^ sin^ ^ 

<1. 



"4cos2^(aj2+j;2) 
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Figure 7. Validity of the quasi-longitudinal 
approximation for the model ionosphere shown 
in figure 3. 

This condition is compatible with the requirement that cos 6 change sign on reflection. Figure 
7 shows the angles 6 for which |/3i^=l and 0.1 computed from the model shown in figure 3. 

It is seen that the quasi-longitudinal approximation becomes doubtful between 75 and 
95 km with our present model ionosphere. At lower altitudes where the collision frequency is 
much greater than the gyrofrequency and at higher altitudes where the electron plasma fre- 
quency is greater than the gyrofrequency, the quasi-longitudianl approximation is useful. 
Since daytime reflection of VLF radio waves usually occurs below 75 km [Bracewell et al., 
1951] and since most of the absorption occurs below 75 km, this analysis is useful for a qualita- 
tive discussion of these effects. 

At altitudes around 85 km the breakdown in the quasi-longitudinal approximation is 
somewhat more serious. For here |/5l^=0.6 for vertical propagation. However, it can be 
seen from figures 4 and 5 that the elements in the coupling matrix C are small in this altitude 
range, so the accumulative effect of the error would not be too great. 

Caution should be used in treating cases with highly oblique propagation, however; the 
quasi-longitudinal approximation is improved by the fact that as the wave propagates beyond 
70 km, its wave normal becomes more vertical. This analysis would break down rather badly 
using a nighttime model ionosphere because there would be a wide range of altitudes at which 

2 

the quantity - .^ ,^ ^ would be considerably greater than unity. However, during times of 

enhanced Z>-region ionization produced by solar flares, auroral, and PCA events, the quasi- 
longitudinal approximation is more generally applicable. Also at higher geomagnetic latitudes 
where the magnetic field is closer to the vertical, the usefulness of quasi-longitudinal approxi- 
mation is improved. 

I thank Miss Estelle Guilbault for her able assistance in carrying out the numerical 
computations. 
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6. Appendix. List of Symbols 

A = matrix [see eq (20)] 

A^ = derivative of the matrix A 

a]^ =[see eq (11)] 

Bq = magnetic induction of the earth's field 

C = coupling matrix =^~^^' 

C^ = coupling matrix [see eq (29)] 

C^ = coupling matrix [see eq (29)] 

Cij = element of C 

c = velocity of light 

c = cosine of zenith angle 

E = electric field strength 

E =E exp ( — ik cos \p3) 

e = electronic charge 

e =base of natural logarithm 

/ = propagation frequency 

k = propagation constant 

m = electron mass 

Ne —electron density 

Q = diagonal matrix of elements kj cos \l/j 

2 = altitude 

a = angle between phase normal and ray direction 

8 =dip angle of magnetic field 

€ = integer function 

6 = angle between direction of propagation and magnetic field 
xp = zenith angle of phase normal 

V = electron collision frequency 

=azimuthal angle of phase normal 

7 = phase angle 

CO =27r/" angular frequency 

co^= ^square of electron plasma angular frequency 

ccm = electron gyrofrequency 

\}c =unit vector in direction of phase normal 

Ift =unit vector in direction of magnetic field 

ly^ =unit vector to specify polarization of electric field (see fig. 1) 

\z' =unit vector to specify polarization of electric field (see fig. 1) 
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